High- Activity Perturbation Expansion for the Hard Square Lattice Gas 
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We study a system of particles with nearest and next-nearest-neighbour exclusion on the square 
lattice (hard squares). This system undergoes a transition from a fluid phase at low density to 
a columnar ordered phase at high density. We develop a systematic high-activity perturbation 
expansion for the free energy per site about a state with perfect columnar order. We show that the 
different terms of the series can be regrouped to get a Mayer-like series for a polydisperse system 
of interacting vertical rods in which the n-th term is of order z~ < -" +1 ^ 2 , where z is the fugacity 
associated with each particle. We sum this series to get the exact expansion to order 1 / z 3 ^ 2 . 

PACS numbers: 75.10.Jm 



I. INTRODUCTION 

Properties of systems with only hard core repulsive in- 
teractions between particles have been of continued inter- 
est in statistical physics. The temperature of the system 
plays no role and the phase transitions in these systems as 
a function of density are purely the result of geometrical 
effects of excluded volume interactions. These have been 
called geometrical phase transitions, the best studied of 
which are the percolation transition |l| and phase tran- 
sitions in assemblies of hard spheres [2[. Many systems 
with different shapes of particles have been studied in the 
literature, for example squares i, rods 0, triangles @ 
and L-shaped molecules Q . The computation of the par- 
tition function of particles with finite exclusion volumes 
on a lattice remains an outstanding problem 0, H| . The 
model of hard hexagons on the triangular lattice and re- 
lated models are the only known exactly soluble systems 
of this kind [t| [lfj. Such models in d-dimensions can 
also be related to directed and undirected lattice animal 
enumeration problems in (d + 1) and (d + 2)-dimensions 
respectively [ill Il2j. 

In many cases, the high density states exhibit crys- 
talline order: the particles preferentially occupy a sublat- 
tice of the full lattice (for example hard hexagons on the 
triangular lattice and the nearest-neighbour-exclusion 
lattice gas on the square lattice). In such cases, it is 
straightforward to develop perturbation expansions in 
inverse powers of fugacity for thermodynamic quanti- 
ties [HJ]. In this paper we study particles with near- 
est and next-nearest-neighbour exclusion on the square 
lattice (hard squares). Equivalently, each particle is a 
2x2 square that occupies 4 elementary plaquettes of the 
square lattice (Fig. [TJ. This system displays "columnar 
order" at high densities where one of the odd or even rows 
(or columns) is preferentially occupied. In this case, the 
standard high-density expansion in inverse powers of the 
fugacity z breaks down. It was realised quite early that 
the leading order correction to the high-activity expan- 
sion is of order 1/y/z fhU l. but a systematic expansion 
has not been developed so far. 

There have been only a few theoretical studies of 
the columnar ordered state so far, and the present un- 



derstanding is not very satisfactory. Correlations in 
the columnar ordered state are hard to capture using 
mean- field like descriptions. For example, all the well- 
known approximation schemes like the mean- field theory 
and cluster variational approximations underestimate the 
value of the critical point z c for the hard squares problem 
by an order of magnitude. 

In this paper we develop a systematic expansion of the 
free energy in the columnar ordered state extending the 
earlier treatment in (l4| . This is a singular perturbation 
series in powers of 1/ \fz. One introduces explicit sym- 
metry breaking by ascribing different fugacities za and 
zb to the particles on even and odd rows. The point 
zb = corresponds to a fully columnar ordered config- 
uration. The perturbation expansion about the ordered 
state in powers of zb is a standard Mayer-like series fl5j |. 
For the case za = zb — z, we show that the terms of 
the series can be regrouped and the resulting series can 
be thought of as a Mayer-like expansion of extended ob- 
jects (vertical rods), but of arbitrary size. The term of 
order l/z~ i~ involves at most n such objects and gets 
contributions from all possible sizes of these objects. We 
evaluate explicitly the contribution of terms correspond- 
ing to two rods, which gives exact results up to order 
1/z 3 / 2 , adding one more term to the earlier known result 
of order 1/z . 



II. THE MODEL 

We consider a system of particles with nearest and 
next-nearest-neighbour exclusion on the square lattice. 
This model has been the subject of many studies. Trans- 
fer matrix techniques indicate the existence of a phase 
transition in this model [l6l - [20j . Variational, density 
functional methods and virial expansions have also been 
used to study this problem fl7l . l2ll . [22[ . These studies in- 
dicate that at high densities the system is not sublatticc 
ordered but exhibits columnar order. However, there is 
as yet no rigorous proof of this. In this context it seems 
worthwhile to develop systematic series expansions for 
the high-density phase, and investigate their convergence 
properties. Recent Monte Carlo evidence suggests that 
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FIG. 1: A typical configuration of hard squares on the lattice 
at (left) low fugacity and (right) high fugacity. The light 
green squares correspond to particles on row A, whereas the 
dark green squares correspond to particles on row B. At high 
densities, the system is in a columnar ordered phase. 
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2 C = 97.5 ± 0.5 [23j, |2J. The nature of the transition 
from the low density disordered state to the columnar 
ordered state at high densities has been the subject of 
several studies. The hard square lattice gas also arises 
as a limiting case of spin models with interactions u p to 
the second nearest neighbour on the square lattice [271 - 
|29| . In the Landau theory paradigm, the four columnar 
ordered states in this system give rise to an XY model 
with four fold anisotropy [30l ] . Such models are known to 
display non-universal behaviour with critical exponents 
that are governed by marginal operators, and naturally 
fall into the class of Ashkin- Teller-Potts models that ex- 
hibit such behaviour j3ll - l33j . However, the critical ex- 
ponents of this model have been hard to pin down from 
Monte Carlo simulations on small system sizes (20j . 

There is some controversy about the nature of the tran- 
sition from fluid to columnar order in this system. Some 
recent Monte Carlo studies suggested an Ising universal- 
ity class [23j], however, this does not seem likely 0, H^]. 
Our own Monte Carlo studies and theoretical arguments 
also imply that the model is not in the Ising universality 
class 12511 . 
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FIG. 2: We define four sublattices 1 through 4 on the square 
lattice. Rows containing the sites of sublattice 1 and 2 (3 
and 4) will be called A (B) rows. Similarly, columns with 
sublattice sites 1 and 4 (2 and 3) will be called C (D) columns. 



the transition from fluid to columnar order in this sys- 
tem is of second order, with exponents close to the two 
dimensional Ising model [23U25I ] . 

We write the grand canonical partition function of the 
hard square lattice gas model as 



^(C) 
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(1) 



where z is the fugacity of the particles, n(C) is the num- 
ber of particles in the configuration C and the sum is over 
all allowed configurations of hard squares on the lattice. 

The Landau free energy per site f(z) is defined as the 
large- N limit of [— log Cl(z)] /N, where N is the total 
number of sites in the lattice. The density of particles at 



a particular fugacity is given by p(z) 



The low-activity Mayer series of this model can be 
computed using standard techniques [l5|. We get (2(| 
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This expansion has a finite radius of convergence de- 
termined by the singularity closest to the origin at z « 
—0.1. The low-density fluid phase has short-ranged corre- 
lations between particles, but the high-density columnar- 
ordered phase has long-range order. Monte Carlo simu- 
lations yield the best estimates of the critical point as 



III. EXPANSION ABOUT THE CRYSTALLINE 
ORDERED STATE 

We consider the square lattice to be made up of four 
sublattices labelled 1 to 4 (Fig. [2]). Rows containing 
the sites of sublattice 1 and 2 are called A-rows, and 
those containing sites of sublattice 3 and 4 are B-iows. 
Similarly, columns with sublattice sites 1 and 4 are called 
C-columns and those with sublattice sites 2 and 3 are D- 
columns. We will specify the position of the square by 
the position of the top left corner of the square. 

At the point z = 00 the problem of computing the par- 
tition function reduces to one of calculating the number 
of perfect tilings of the lattice by hard squares. On a lat- 
tice of size LXL (where L is assumed even), with open 
boundary conditions, clearly 



a. 



,{z -> 00) - z l2 / 4 . 



(3) 



For a lattice with periodic boundary conditions in 
both x- and y- directions, it is easy to see that we have 

Aperiodic — 4(2 ^ -0^ ^ ■ 

For z finite, there is a finite density of plaquettes not 
covered by squares in the system. Suppose we start with 
a fully packed configuration of squares, with all sites on 
sublattice 1 and introduce a single vacancy. By sliding 
squares in the corresponding row, we can break the va- 
cancy into two half-vacancies with in-between squares on 
the 2-sublattice. These two vacancies can be moved ar- 
bitrarily far apart. Thus there are 0(L 2 ) configurations 
of half- vacancies on each row, and there are L/2 possible 
rows. There is an equal number of configurations with 
half- vacancies along columns. This leads to 
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This phenomena, which may be called deconfinement 
of half- vacancies, is the reason behind the failure of the 
standard Mayer expansion technique, in which the order 
1/z term has a coefficient of the order of the volume of the 
system. Because the standard cumulant expansion fails, 
this series must be treated using different techniques. 

A simple qualitative picture of the system with a small 
density of defects introduced on the background with per- 
fect crystalline order is this: The pairs of half-vacancies 
may be pictured as joined by straight rods, which may 
be horizontal or vertical. The state with oricntational 
(i.e. nematic) ordering of these rods corresponds to a 
state of the hard squares with columnar order. In fact, 
there is an exactly solved model of hard rods of variable 
length due to Ioffc et al. (34J, in which the activity of 
a rod is independent of its length, which shows a phase 
transition at a finite value of the activity of the rods. 
This phase transition is in the Ising universality class. 
However, we note that the precise weights in Ioffe et al.'s 
model arc different from the hard squares case, and, in 
fact, the transitions are not expected to be in the same 
universality class. 



IV. EXPANSION ABOUT THE COLUMNAR 
ORDERED STATE 

In the columnar ordered state the system preferentially 
occupies one of the rows (A or B) or columns (C or D). 
In the A-ordered phase we have (pi = p-i) > (p3 = Pi), 
where pi denotes the density of particles corresponding 
to the ith. sublattice. The B, C, and D phases are defined 
similarly. 

To quantify the nature of ordering in this system we 
define the following order parameters. The row order 
parameter of the system is defined to be 



R = A[(p 1 + p 2 )-(p 3 + p4)}, (5) 
and the column order parameter is 

Oc=4[(pi+p4)-(p2 + Ps)]- (6) 

Equivalently, we can also define a single Z4 complex 
order parameter 



iV2[(p t - p 3 ) + i(p 2 - Pi )}. 



(7) 



The factor 4\/2 has been introduced to make the 
maximum value of the order parameter 1. The phase 
of the complex order parameter Oz 4 takes the values 
7r/4, — 37r/4, -7r/4 and 37r/4 in the A, B, C, and D phases 
respectively. 



We now develop a high-activity perturbation scries 
about the row-ordered state. We effectively integrate out 
the horizontal rods, and this generates a longer-ranged ef- 
fective interaction between the vertical rods, in addition 
to the excluded volume interactions. 

As in [3] j we associate different fugacities to the parti- 
cles on the A and _B-rows. The fugacities of the particles 
on row A is za and on row B is zb- FI(za, zb) is the par- 
tition function of the system with this explicit symmetry 
breaking between the two rows. H(za,0) corresponds to 
the fully columnar ordered state. We then expand the 
partition function of the system about a state with per- 
fect columnar order in the A-phase. 

We write the partition function of the system as a for- 
mal expansion in terms of the fugacities of the particles 
on the -B-rows (defects) and the corresponding partition 
functions of the A-rows. The partition function expan- 
sion about the columnar ordered state is 



^^ = 1 + z b W 1 (za)+ Z -^W 2 (z a ) + .... (8) 
u{za,0) 2! 

Now, taking the logarithm of Eq. © , we arrive at the 
cumulant expansion 



zb 2 

\ogVL{z A ,z B ) = ko(za) + zbki(za) + -^-K2(za) + - 



N 



(9) 

The calculation of W n (zA) involves fixing a configu- 
ration B n of n particles on the S-sites. The weight of 
this configuration is defined to be zgProb(B„), where 
Prob(S„) is the probability that in the reference sys- 
tem with only A-particles, no site excluded by these B- 
particles will be occupied. We then sum these weights 
over all configurations B n . The negative of the logarithm 
of Prob(i3 n ) is defined to be the effective interaction en- 
ergy between the B-particles when the A-particles are 
integrated out. 

It is straightforward to evaluate the first few terms in 
this expansion. When there are no -B-particles in the 
lattice, the partition function of the system breaks up 
into a product of 1-d partition functions of particles on 
the A- rows (particles on different A-rows do not interact). 
These A-particles thus behave as a 1-d lattice gas with 
nearest neighbour exclusion. Thus we have 



n(* A ,o) = [n ldtL ( ZA )] L/2 ■ 



(10) 



where ^ici,l(za) is the 1-d partition function of particles 
with nearest neighbour exclusion on a periodic ring of 
length L. This is easily seen to be 



n ld , L (z A ) = \+ L + \- L , (11) 

where A± are the eigenvalues of the 2X2 transfer matrix 
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z A 
with 



A± = 



1 ± VI + 4za 



(12) 



In the limit of large L, this gives us 



k-o(za) = 2 lo 8' A +- 



(13) 



At the next order, W\{z A ) involves fixing the position 
of a single particle on a -B-site, say at {x,y). Then the 
sum over A-particlc configurations is restricted to those 
in which the sites (x ± 1, y ± 1) and (x, y ± 1) are not 
occupied. The partition function of a single row with 
this constraint is the partition function of a system of 
particles on an open chain of length L — 3, and is easily 
calculated. Let /ooo be the probability that the sites in 
a randomly picked interval of length three on an A-row 
are all empty, in a 1-dimensional lattice gas with nearest 
neighbor exclusion at activity z A . Then, 



N , 
Wi(za) = y (/ooo) 2 . 



(14) 



The expression for /ooo is easy to derive using the prop- 
erties of the 1-d nearest-neighbor exclusion lattice gas. 
We have /ooo = j^foio, where /oio is the probability that 
a randomly picked site in the gas will have the occupa- 
tion numbers 010 at the three consecutive sites. Clearly, 
/oio = Pi<i(za), the density of the gas. Using the Eq. 
(fT2"|. we get 



P\d{z A ) = -- 



2 2 VI + 4za 



(15) 



Therefore we obtain 



ki(za) = ^ (/ooo) 2 = ^ 



Pld(zA) 



ZA 



(16) 



Therefore the leading contribution to the cumulant ex- 
pansion from the single particle term is of order O {jf^j ■ 

The calculation of higher order terms W n (zA) for n > 2 
is similar. However, working order by order in n is not 
very effective. In the series given in Eq. ©, the series 
is in powers of zb, with coefficients that are functions 
of za- Eventually, We would like to put za = z B = z, 
and expand the series in inverse powers of 1/z. Unfor- 
tunately, the leading behaviour of W n (zA) for large za 
is only z A ^ n+1 \ Then, arbitrarily large orders in n are 
required even to get the correct result to order 1/z [lij . 



This may be seen as follows: Consider the configura- 
tion of n £?-particles vertically above each other. It is 
easy to see that the probability that such a configura- 
tion would be allowed in the unperturbed ensemble is 

[/ooo]" +1 , an d hence is O ^z^" +1 ^. Hence the contri- 
bution of this term to the perturbation expansion is of 
order zb" z^ n+1 K For zb = za = z this is 0(1/ z) for all 



It is easy to check that all configurations except these 
vertical rod-like configurations of £?-particles do not con- 
tribute to order 1/z. Hence, if we want to sum over all 
terms to order 1/z, we group these configurations to- 
gether, and identify them as a single vertical rod. A gen- 
eral configuration of B-particles would then be a group 
of non-overlapping vertical rods. In any configuration of 
-B-particles, we define a rod containing a given occupied 
B-site f to be the set of all the consecutively occupied B- 
sites in the same column reachable from r using vertical 
steps of length 2. Clearly, any configuration of B-sites 
has a unique description as a set of non-overlapping ver- 
tical rods. To avoid over-counting, two rods cannot sit 
directly on top of each other. 



V. HIGH-ACTIVITY EXPANSION IN TERMS 
OF RODS 

We now develop an expansion in terms of the number 
of rod defects. We attach an additional activity factor e 
to each rod, and now rewrite the summation in Eq. (|5]) as 
a sum over configurations involving different numbers of 
rods 



Q(z A ,z B ) 

n(z A ,o) 



= 1 + eRi(z A , z B ) + e 2 R 2 (z A , z B ) + ■ ■ ■ , (17) 



where R u (za, z B ) denotes the contribution of the con- 
figurations with exactly n rods to the partition function 
expansion. Now taking the logarithm, we arrive at the 
free energy expansion 



T(z a ,Zb) 



1 

"N 



log (Cl(z A ,z B )) ■ 



= F (z A , 0) + cF^Za, z B ) + e 2 F 2 (z A , z B ) + . . . , (18) 

where F{z a ,zb) denotes the free energy per site of the 
hard square lattice gas and 



F\(z A ,z B ) = -—R^za^Zb), 

VI \ 1 I U I \ Rl(zA,Z B f 

F 2 (za,z b ) = -— I R2(za,zb) 2 



; (19) 

We evaluate the free energy expansion formally in pow- 
ers of e. At the end of the calculation, we set e — 1. Fq 
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denotes the contribution from the term when there are occupied by both rods (The weight of the configuration 
no rods in the system. From Eq. (fT5|) we have w(fi, h;?2, h) is zero, but w(fi, h)w(r2, 1 2) is finite)(Fig. 

©. 

Fq = — k o( z a)- (20) We deal with these three terms separately. We have 



F2 — t\ + t2 + £3 



(25) 



FIG. 3: The configurations contributing to a single rod term 

Now consider configurations with exactly one rod. 
The weight of a rod with length n is easily seen to be 
/ooo" +1 z B n ■ There are Nb = N/2 possible positions for 
each size n of the rod. Summing over all possible values 
of n, we get 



N 



which yields, 
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It is convenient to define a parameter 7 
Then we have 
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Wc next evaluate the two rod term, with 



(21) 

(22) 
z bJooo- 

(23) 
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fi,h f 2 ,h 



w(ri,h;r2,l 2 ) 



w(ri,h)w(f f 2,h) 



(24) 

where w(ri, l\\ r*2, h) denotes the weight of the two-rod 
configuration with the centre of the first rod of length l\ 
at position f\ and the centre of the second rod of length 
I2 at position r 2 , in the partition function expansion. 
w(fi, h) denotes the weight of a single rod configuration. 

In this summation, the contribution from the two-rod 
terms that have zero interaction are exactly cancelled 
in the corresponding cumulant expansion. The configu- 
rations with non-zero interaction between the rods are 
ones where at least one row is touched by both the rods. 
These configurations can be classified as follows: 

1) Adjacent rods (configurations in which only one row 
is touched by the ends of two rods from opposite sides 

(Fig. a), 

2) Two rods which share a non-zero interval in the 
Y-direction. In this case there are a finite number of 
rows occupied by both rods (Fig. [5]), 

3) Overlapping rods, where a finite area of the lattice is 



1. Adjacent rods 



h = 2 i.. 



FIG. 4: The configurations contributing to the adjacent rods 
term 

In the calculation of the term ti , we deal with all con- 
figurations in which two rods touch an A-row from dif- 
ferent sides. Then for this row, we need to calculate the 
correlation function in the reference problem that two 
triplets of three consecutive sites at a given distance are 
both empty. This is easily expressed in terms of the cor- 
relation function of occupied sites in the reference sys- 
tem. Let the pair correlation function at a separation 
A between the centres of the triplets in the f-d nearest 
neighbour exclusion gas be G(A). Then it is easily seen 
that 



G(A) = (jooo) forA = 0, 

= (/ooo) ~r~ for |A| = f , 

= (/ooo) 2 (l - « |AM ) for |A| > 2, 



where 



a 



-A 



1/2 



o 



3/2 



(26) 



(27) 



One of the rods extends upwards with length n a and 
the other rod extends downwards with a length rib (Fig. 
2J, the factor of 2 accounts for the symmetry related dia- 
grams. Summing over all configurations and subtracting 
the disconnected part, we have 
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(28) 
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This scries can be easily summed. We have 



7 2 /ooo / 1 /( 



ooo 



(1-7)2 VA+ I- a 



(29) 



2. Rods which share a non-zero interval in the Y-direction 



FIG. 5: The configurations contributing to the term involving 
rods which share a finite Y-interval. The factors multiplying 
the diagrams account for the symmetry related configurations. 

In the evaluation of the term t%, we have two rods 
with n i?-particles sharing a common interval in the Y- 
direction. Now, one of the rods can have a part extending 
above this section (Fig. [5]), this can happen in 4 ways. 
Also, the two rods can both have a finite extension above 
this section, this can occur in two distinct ways as follows: 
one of the rods extends both above and below the other, 
or one rod extends above and the other below. Both these 
cases yield the same weights. The extending sections 
above and below the overlap are of lengths n a and rib- 
We have 

1 / oo oo oo \ 

t2 = ~- l + 4^ 7 "»+4^ J2r a+n ") 

\ ri a — 1 n a — l n^ — l / 

oo oo 

£ E7 2 "°(/ooo) 



, n a = l A>2 



i_ a .Ai-iy* +1 -i 



(30) 



The sum within the first brackets can be performed easily 
to give fjz^J • The summation in the second brackets 
can be evaluated and we have 
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7 2 (/ooo) 2 



, A>2 
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(/ooo) 7 
1 - 7 2 



(31) 



We can rewrite this in a better form as 
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Now, since a is negative, the terms in the summand 
series are oscillatory. We thus split the terms into even 
and odd powers of a. We write 



ta = 
where 

Seven ' 

and 

Sodd = 
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(33) 



(34) 



(35) 



For large z, a ~ — l + z -1 ^ 2 , and the summand in S dd 
and S even is a slowly varying function of n. To evaluate 
these expressions we can approximate the sum over the 
discrete values by an integral over the appropriate range 
of n. The error is only of order 1/ z. We get 
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5. Overlapping Rods 
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FIG. 6: The configurations contributing to the overlapping 
rods term 



In the calculation of the overlapping rods term, we deal 
with all possible configurations in which two rods which 
partially, or fully overlap. In this case the two rods can 
be at a distance A = or A = ±1 away from each other. 
The length of the overlap is n and the parts extending 
above and below are of lengths n a and rib respectively. 
There is also a contribution from the forbidden configu- 
ration in which two rods sit directly on top of each other. 
We have 
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n a — 1 / n a —l — 1 



+rib 



which can be easily simplified to give 



ta = - 



3 /l + 7\ 2 //ooo 2 7 2 A , 1 //ooo7 X 2 



4 V 1 ~ 7 / V 1 - 7 / 2 V 1 - 7 
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A. High- Activity Expansion for the Free Energy 



Setting e = 1 in Eq. (|18|) we obtain a high-activity 
expansion for the free energy of the system. The term F\ , 
corresponding to single rods, has a leading contribution 
of order 1/z. The term F2, corresponding to two-rod 
configurations has a leading contribution of order 1/z 3 / 2 . 
The order 1 / z 3 / 2 contribution to the two- rod term comes 
from two sources, a) when two adjacent rods separated 
by a distance 1 touch the same row from opposite sides 
and b) a summation over the distance between two rods 
which share a finite interval in the Y-direction. The sum 
over the distance in this case yields a factor of order y/z 
because the correlation length between the rods is £ = 
log a ~ y/z. The term of order 1/z 2 gets contributions 
from configurations involving three rods, in addition to 
the above single rod and two rod terms. For n rods, 
the summation over the distance between rods yields a 
leading order contribution of at most z~ ~ . Hence in the 
evaluation of the term of order z~?~ , only terms involving 
n rods need to be considered. 

Using our expressions for Fq, F± and F%, we can gener- 
ate the exact series expansion for the free energy and the 
density of the hard square lattice gas up to order 1/z 3 / 2 . 
We have 



-/(*) = 

4 4Z 1 / 2 4z 



.3/2 



and 

1 1 1 (flogf 9 ^ ' ir 
4z 



O 



(40) 



P{ ' 4 8zV2 A , 



,3/2 



o\- n 



(41) 



Higher order terms can in principle be evaluated in a 
similar manner. The term of order 1/z 2 gets contribu- 
tions from configurations with up to three rods. This in- 
volves many more diagrams with more complicated sum- 
mations, and will not be attempted here. However the 



convergence of this series for z > 150 to the actual val- 
ues seems quite good. For example at z = 200 and 150 
Eq. (gTJ gives p = 0.2396 and 0.2377, to be compared 
with results from Monte Carlo simulations which yield 
p = 0.2395 and 0.2374 respectively. 



B. Order Parameter Expansion 

In this section we compute the high-activity expansion 
for the row order parameter which is defined in Eq. ([3]). 
We have 



O r (z a ,z b ) = 4 



2,1 



dz A 



ZB 



dz B 



(-F(z A ,z B )) ■ 



(42) 

Using the terms in the free energy expansion above 
we can generate the exact order parameter expansion to 
order 1/z 3 / 2 . We have 



Or(z, z) 

1 



1 



2^/2 



395 „„, /9\\ 1 , n \ 



(43) 



Using the terms up to order 1/z 3 / 2 in the order pa- 
rameter expansion, we estimate the critical point of the 
system by solving for the point On = 0. We obtain an 
estimate of the critical point z c = 14.86, which is signif- 
icantly better than the estimate z c = 6.25 obtained by 
using just the first three terms of the series. 



VI. DISCUSSION 

In this paper we have developed a perturbation ex- 
pansion about the columnar ordered state for the hard 
square lattice gas. We identified the basic objects of ex- 
citations about the ordered state, namely vertical rods. 
We showed that only configurations with at most n rods 
contribute up to 0(z~( n+1 ^ 2 ) in the expansion in inverse 
powers of z" 1 / 2 , and we explicitly summed the contribu- 
tion from the terms containing two rods that provides 
exact results up to order 1/z 3 / 2 . It is possible to extend 
the series expansion developed in this paper to other sys- 
tems that display columnar ordered behaviour at high 
densities. In particular for the k x k hard square lat- 
tice gas a similar procedure can be used to generate the 
high-activity expansion with terms of powers of z _1 / fc 
appearing. Other systems displaying columnar ordering 
at high densities such as hard rectangles on the square 
lattice are also amenable to similar treatments. 

The method can be extended to three dimensional sys- 
tems that exhibit columnar order. Consider a system of 
2x2x2 hard cubes on the cubic lattice. Consider a state 
with perfect columnar order, with columns along the z- 
direction. Then its cross-section in the xy-plane would 
be a fully packed configuration of 2 x 2 hard squares. 
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Then, as discussed above, there are approximately 2 L / 2 
different fully-packed hard squares configurations. We 
can develop a perturbation theory in inverse powers of 
y/z, about any of these states. However, it is easily seen 
that the first correction term (of order 1/z) is different for 
different states, and is largest when all the columns are 
arranged in a periodic supcrlatticc with square symme- 
try, and spacing 2. These states then outweigh all other 
states, for nonzero 1/z, in the thermodynamic limit. This 
is an example of "order-by-disorder" [35[ in this system. 
The extended objects that contribute to order 1/z in the 
za — zb = z series in this case turn out to be rigid rods 
along the x- or y- directions. The interactions between 
the rods are similar to the two dimensional case, but the 



summations are harder to do in closed form. This seems 
to be an interesting direction for future studies. Another 
problem of interest is to establish lower bounds on the 
radius of convergence of these expansions. 
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